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ABSTRACT 


The Wavenumber Technique (WI) is a relatively new method 
of underwater sound transmission analysis. One aspect, 
source, depth detenai nat ion, 1s studied to evaluate its 
validity and test environmental and acoustic sensitivity. 
The horizontal wavenumber spectrum is analyzed to determine 
null spacings in wavenumber space, which indicates source 
depth by the Lloyd's Mirror interference affect. Comparison 
of this theory with cases of an isosp2ed sound profile, 
Pel ly S@hbserp2=ng bottod, and flat totally-reflecting surface 
shows excelien* agreement for several parametric variations. 
Cases with realistic sound speed profiles and partiaily 
absorbing bottoms genérally agree with theory, but a dis- 
tinct bles is observed. Source cepth determination curves, 
which relate tne scaled wavenumber spectral intensity null 
spacing to the source depth, are presented for comparison 
W2Eh theory and sensitivity analysis. An example is given 


for suggested application of source depth determination. 
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The Wavenumber Technique (WT) is anew method for the 
solution of several problems related t9 the analysis of 
underwater sound transmission (Lauer,1979). Its basis is the 
analysis of acoustic propagation in the wavenumber domain, 
which has been described by DiNapoli (1971) and was an 
antexrmediate step in the computation of transmission loss 
versus range in the Fast 252d Program (FFP). DiNapoli 
(1977) later used the wavenumb2r domain to study the impe- 
dance of the ocean-bottom intzee face. Lauer (1979) first 
heemegucec the WI and its applicatwons t2 voassive localiza- 
erom aid tracking and multipath decomposition. He daveloped 
a rcange prediction curve based upon peak wavenumbers in the 
transform spectrum fcr successive range increments from an 
entire pressure range field. 

These investigations presented highly beneficial and 
promising results from the wavenumber depiction, which indi- 
Gast-oueM@ecmecc: <20na@il frssearch is appropriate. One particu- 


the WR - determination of source depth - will 


cr 
O 
ty 


a 


rey) 


mr aspec 
be investigated here in order to evaiuste its vaiidity in 
Go@perison with «the Lloyd's Mirror eéisct. Ets ela. on- 
Sae0s 25 oceanic and acoustics variations will also be 


= 


examined. 
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The wavenumber technique as first described by Lauer 
(1979) was based upon input from the FFP, but it has been 
adapted for use with the Parabolic Equation (PE) method of 
propagation loss determination. The PE-based WT is currently 
being investigated by Lauer and others at the Naval Ocean 
Research and Development Activity (NORDA) and was the nodel 
used in this research a* th2 Naval Postgraduate School 
(NPS). Descriptions of the PE and WT will be presented to 


display their sensitivities to variations of several model 


oO 


14 


atameters. Also, ¢xfected relationships of the WT to other 
current propagation prediction models will be discussed «+o 


consider naow resuits might vary based upon different inputs. 


Pate eteular attention will be given to algorithm differences, 


de¢qrees of approximation, and =réatment of «he ocean botton. 
The E-pased WT will then be exercised for several 
parameter variations to evaluate qualitativaly the consis- 


eewrepect wRe technique. This study will also demonstrate 


mae Seeponse of the PE-based WI to geometric and actual 
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es 


paiee 86 COG Lticns. Pred:¢z292n curves fo2 source depth 
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A. GONERAL DESCRIPTICN 

The WI is a natural alternative for signal depiction 
because it prcvides information on the directions of energy 
arrivals. I+ allows reasonably direct physical interpreta- 
tion of ocean acoustic processes and their relationships to 
environmental conditicns (Lauer, 1979). Lauer further states 
that «he wavenumber spectrum plot is clearer and more 
informative than typical curves of propagation loss versus 
range. The wavenumber K relates the angular frequency and 


he sound speed by 





= ae ae 
— J 
las Qn 
VW ~ 
pede 2st mes hepazontal and vertical componen«s 


The reference wavenumber is determined by the sound speed 
AL aAigun 
K = 0 
QO =~ 
C 





The WI wiil first be described without reference +o any 
particular propagation model t> present the general nature 
Of the algoritha. Fig. Tis the flow diagram for the WT. 
Meee Welt ahd right sides of this flow diagram correspond “<o 
the intensity and wavenumber axes, respectively, of the 
Seeeccem plot’ that will result from the Wf. This figure and 
the following discussion were developed from Pamir (1979 and 
1982a) and augmented by several telephone conversations with 
Re. Lauer, T. Lawrence, and R. Evans of NORDA between July 
and December, 1982. 

The WP c@iculates spectral intensity from the real and 


imaginary parts of.the acoustic pressure 
Lut 


SS Ot) a ; 


The complex psi versus range field at 323 specified depth is 


the propagation model product which deteraines the pressure. 


tly 


zed 


hh 


tt 


The compiex pressure field can be nodi Cie gertects such 
és attenuation, depending upon whether liealized or realis- 


= 
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fic comditions ars teing investigated. A Fast Fou 
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S-amc"Grm {PFT) is then applied to yield transformed complex 


pressure. The spectre ‘ntensity is determined a+ #ach 
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wavenumber increment from the transformed complex pressure 


field by 


ce a 2 
I Re im, : 


This intensity field is scaled arbitrarily for satisfactory 
display of the relative maxina and minima to define the 


wmilis. 
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Additional output parameters from the PE which aré 
Stequired for the wT include th2 range stao, the number of 


Tange points, frequency, Minimum sound speed and field 


og 


. These are the seiected variables noted in figure 1. 


fu 


S 


'O 


- 
- 


3 


he range step is selected to provide sufficient resolution 


=o determine =he wavenumber null spacing on the spectrum 
DLOtws. The initial radial component of the wavenumber for 


ct 


me opee-cum -5 Cadculated from DiNapoli (1971) by 
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AS will be seen in section II-B, if a "scaled" wavenun- 
ber beta is used for the horizontal axis, analysis is sin- 
plified since the nulls in the spectra shouid be more evenly 
spaced (Lauer,1979). The physical meaning of beta will be 
described in the next section. The resulting scaled wave- 
number spectra are then plotted. The wavenumber null Spacing 
tor each incremental source depth is determined from each 
scaled spectrum. The securce depth determination curve is 
generated by plotting these null spacings as a function of 


source depth, as suggested by Lauer (1979). The nuil spac- 


t 


ng cf an analyzed received signal could be compared with 
Ehis GewenWination curve to infer the source depth of the 
signal. This application will be discussed later. The 
Source depth determination curve will b2 specific for the 


Peomeetc and oceanic descriptions of the mediun. 


B. Peis l@@e DESCRIPTION 

The precéding description was intended to provide the 
ecverall concept of «he WT whersas whazt follows will explain 
Pie wimwetcat basis for che technique. The underlying princi- 
pia cf the WT as descriped by Lauer (19382a) is the Litova"s 
reborn ittec., Which giyes the acoustic field by 


a = PC) aleee F 


PAS 





where 


ECR) = exp {ik / R2 + ree) - exp tik / R2 + Ca eM 


/ 2 e 2 / 2 2 
R< + CZ. Zo) R*+ + (Z. + Zo) 


for a point source ina semi-infinite nedium of constant 


sound speed with a pressure release surfaces. This is the 


Lloyd's Mirror field in range space, whereas the Lloyd's 


Micrer field in wavenumber space, Jpop icaple to ~eheas 


research, is given by Lauer (1982) as 


and 


ox ee sin(8Z,) exp(i8Z,) 5 


E(R) and P(K) are related by the Bessel transform Da aan 


ECR) = ee F(K) J (KR) K dK 


F(K) 


it 


ECR Cie ioe dR 
9 O 


Se) | j—4 
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For this tesearch the Bessel function was approximated as a 
trigonometric function so that the FFT could be utilized. 
Maas 2S cOms@dered acceptable since tha source and receiver 
are greatly separated, according to Lauer (1982b). Inspec- 
ten Of Fi K) reveals that the nulls of the spectrum will be 


equally spaced at intervals of 


AB 


i a 


This equal spacing in beta is the key element of the apvli- 
cation or the WT to source depth determination (Lauer,1979). 
Only the direct and surface-reflectei waves will inter- 


West tec LOll vO eae 


od 


Mee @ea2 Lrecez yer location in 4 Simple il 


his interec- 


r3 


Peoga™s Witrror effect with a s@ooth surf2ce. 


TLOR Will be destructive or constructive as, a JEaculsecoe 
g 


tment. selatave phases (Kinsler 2t al ,1982) and will produce 
mulls or peaks, respectively, in the wavenumber spsctrun. 
meee, 2 220 2Cces Ere Licyd*s Mirtor effect yeometry which will 


be used +o study test case outputs for null svacing and sub- 
Sequent source depth determination. The surface-reflected 
(SR) Mave clave: distance wiil increase at a faster rate 


shan the direct (D) wave distance as th2 source is moved 
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peomee 22 Licyd's Mirror sffect Geometry. 


PAS) 





G@eeper om the water cclumn (i) while the receiver depth (h) 
2S stationary. This leads to different arrival patterns and 
nulls at more wavenumbers. Thus the general pattern of 
decreased spacing in beta with increasing source iépth is 
observed. 

Bottom interacting waves will be neglected for idealized 
cases. This wiil be accomplished by the use of a fully 
absorbing bet*om. These waves will be important, however, 


in shallow water or at qreater ranges. A few test cases w 


j-- 


aed 
have realistic bottoms. Various propagation models handle 
she bottom in different ways, so the inzlusion of a botton 
will be discussed later. Cases in shallow water (less than 
5000 feet) or at extended ranges (greater than 100 milé¢s) 
will not be considered here. 


Lauez (1979) observed that th? wavenunbeé> spéec=rum shows 


AY 


series of modulating envelopes which enclose spikes that 


he eigenvalues associated with the nsrmal modes. Con 


ray) 
rf 
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arison of the wavenumber spectrum with typical transmission 


loss curves reveals that the wavenumber depiction is more 


COmemeerc sy, structured. This leads +5 a better physical 
understanding of signai transmission because cf the rela- 
tionship between horizontal wavenumber aid tay <heory. de 


gives this as 
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The horizontai wavenumber decreases as theta increases. Thus 
the wavenumber axis of the spectrum plot will include wave- 
numbers associated with smaller elevation (or depression) 
angles in increasing crder of wavenumber (Lauer,1979). se 
is then expected, for the case of a neglsct2d bottom, that 
spectrum ampiitude peaks could be observed that would corr?2- 


spond to the direct and surface reflected waves. 


oF MAVENUBBER TECHNIQUE APPLICATION 

Fawer (1979) recommends a receiver that is a single 
Omnidirectional hydrofhone and 2 scurce-generated continuous 
wave (CW) received signal for source depth determination. 
He further states that in actual use th2 single hydrephone 


could be replaced by an array to improve the signal-to-noise 


c(t 


ratio and to cbtain bearing and bearing rate information. 


An example of a pessible application of the WT, scurce 


t 


Gempth deter@ination, is presented im Fig. 3. A propagation 
model could be run for a series cf incremental source depths 


fo #22 1g te null spacing relationship, as was shown in Fig. 


hen 


Ts Miatecucve would be speciisze to the prorerties of the 


(p 


ocean and values for parameters such as iIrequency and 
receiver depth. A received signal as shown on the tight side 


Pees cOulo then be analyzed. I*s null spacing could 5é 
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WT APPLICATION: SOURCE DEPTH DETERMINATION 
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JT Appiication for Source Depth Determination. 
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Gompapei to the wmodel-produced curve to infer the depth of 


the source. 


D. TSCHNIQUE LIMITATIONS 

The success of the WT is related directly to the accu- 
racy of the predicted pressure field prodjiuced by the propa- 
gation model. There are several possible differences 
betweer propagation models which might affect the pressure 


fiscld. Some of the more significant include approximations 


ri 


or the wave equation, bottom affects, initial source func- 
tions and input limitations. These will be compared for sev- 
eral nodels to determine relevant restrictions. The pressure 
field will also respond to changes in th2 ocean. A series 


Ge Gieeerent scenarios will be run to provide an initial 


Sszgmmere Of the sensitivities to some of these factors. 


}- 


Any additional BAwltati@ns will be rslated to computer 
PCOCeSsincg time and storaqe requirements. These factors will 
not bre addressed directly in this researsth because the pro- 
cessing following the propagation model is minimal. There is 


aiso no requiztement tec store the complex pressure fielas cr 


forms once the nuli spacing has been determined. 


wmes - “fan 
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eas sie WT computer weiiization is directiy proportional to 


she rime and storage required for the propagation models. 
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III. PROPAGATION MODELS 


A. INTRODUCTION 

Modeling the propagation of sound in the ocean is con- 
plicated by the variability of the ocean, the great range of 
frequencies of interest, andthe many applications of sound 
propagation (DiNapoli and Deaven port, 1979}. It is not sur- 
prising that no single model is presently capable of 
addressing ali of these variations; there are many semi-re- 
Strictive methods. In general all models consider the sound 
speed ~9 be a function of the spatial coordinates and inde- 


pendent of «ime (DiNapoli and Deavenport, 1979). Only the 


t ty 


requency, water column geometry and oteanic description 
Will affect propagaticn. 

Ali models, regardiess of application, can be divided 
into two classes: Tange independent or range dependent. 
mange independent models assume that the oc2an is cylindri- 
cally symmetrical, the speed Sf sound is a function only of 
Ge¢pth and ali boundaries are planar perpendicular *o the 
Gep=h awas (DiNapoli and Deavenporce, 1979). Range dependent 


models can provide better approximations of real conditions 


a 








by allowing the sound speed to be a function of two or three 
spatial variables. Teimalisé do not require planar bounda- 
ries perpendicular tcthe depth axis, thus cylindrically 
symmetric bottom topcegraphy can be included (DiNapoli and 
Deavenport,1979). A range-dependent mod2l1, the split-step 
Fast Fourier Transform PE (SSFFT) was the numerical model 
used in this research to generate the pressure fields for 
*he wavenumber analyses. This PE will be qualitatively con- 
pared to two other range-~deperndent aigorithms, a range-modi- 
fied Normal Mode (NM) pene and a Finite Difference (FD) 
model. Zt will alse de compared to @ range~indépendent 


Bezthcd, tne Fas* Field Progran (FFP) model. The range-de- 


fa» 


Pemcen®t Capability cf the SSFFT was not used in this 
Geeqgares Since this was an initial investigation. Its com- 


parison with both range depenient and independent models 


WD 


may, however, indicate possible model efrects on the WT 
Patch Couic be studied later. 

In general ail models attempt to solve the acoustic wave 
=quation. This equation is valid for ssund propagation in 
E%uid media weet h specifi ed boundary Gond2 ti one 
@iz2ok,1975,p.114). Tt becomes the redaced elliptic Hein- 


hoiltz equation in spatial ccordinates under the assumption 


o= narmonic time dependence. 
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B. PARABOLIC EQUATION 
1. Introduction 
The PE algorithm used in this research was isveloped 
py Breck (1978) . Me Has Cc Ledeure sf “his PE is the 
~eplacement of the recuced elliptic Helmholtz equation by a 


paeemoltc partial differentiai equation that utilizes th 


(D 


D 


Daypppert-Hardir split-step Fourisr algorithm (SSFFT) for th 


Newm@rsicail integration (Brock,1978). The 9 


| 


iginai somputsr 
code as described by Brock has been nodified to incorporate 
Zhe source ard a*tenuation descriptions of Tatro (1977) and 
to Mreduece an ocean bottom by Stieglitz at a) ogo) Phe s 
Meeee me WerezOorn Of the PE that is currently available az NPS 
and was the code (PEMODEL) use2i for this research. {It must 


be emphasized that some difficulty was encountered in trying 


RO eme@rteci th> correct d2scription of the NPS versi 


O 


n 
because of inadequate and erroneous comment. card documenta- 


Rien in thé source ccde and incomplete manuscript descrip- 


: ae . ; f = nee 

Meas sec 20n well discuss the SSFFT, the chad? ive 
Beeck @lgetithm, modifications to that coi2, and subsequent 
EVatUerticns of the SSrFT. Particulac atczention will be 


Beer BO ene=e dspects which ars expected tc aftect the WT. 


5)6. 








2. Parabolic Approximation 








The parabolic approximation to the wave equation was 
first introduced by Leontovich and Fock (1946) as a solution 
for the propagation of electromagnetic waves along the sur- 
Beace of the earth. £ts first dmeroduction to wunderweter 
acoustics was by Hardin and Tappert (1973) who used the 
Split —-scep “cehnigue in conjunction with the FFT to solve 
Uaes Be eticigar apprcximation tothe wave equation. The 
result was the SSFFT. Tappert (1977) related the use of the 
parabolic approximaticn to sound channel propagation ina 
Wavegusae that is thin vertically and elongated horizontally 
to the range cf the first convergence zone or farther. 

Demwuetenge propagatien, with which the SSFFT is 
Dasicaliy concerned, will usually be for iow frequencies 
because volume absorption increases strongly with frequency 
(fappert ,1977). The maximum elevation (and depression) 
angles of propagation for lona range transmission must be 
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SQhii 40 satisfy the FE (Tappart,1977). The derivaticn of 


ct 


h 
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SSPPT W2li be discussed iater an a simplified version as 


= 


(iD 


mead By Brock (13978). 


‘oO 


S 


M 


h- appeses twat the SSFFT is intended for usé in 


rey) 


Situations where bottom interasticn and surface scatteri 
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would be negligible. Thus the use of the SSFFT to examine 
chee Lioeyd's Mirror effect in the real world may be inappro- 
priate if the ocean surface is sufficiently rough; addi- 
tional errors may occur in situations where bottom effects 
eee Gignizicant. 
3. Brock Alcorithn 
@. Description 
m2 SSEPT 1s valiil for acoustic pressure ina 
medium of constant density with a monofrequency source and 
cylindrical symmetry about the iepth axis (Brock,1978). The 


acoustic pressure can be written as 


= i 
ae) = eka) CK RE, 


Whezte€ tne zero-order Hankel function of the first kind 
selates the acousti pressur2 +¢+0O an outward propagating 
Seer cel Weve envelope function (Brock,1978). This is 
valid because of =the assumption that at low frequencies all 


icant energy will propagate approxinately horizontally 


bth 
po 


eiqn3 


away from the source (Tappert,1977). The asymptotic form of 


the deakei function is 





if the receiver is many wavelengths from the source. This 
leads to (Brock, 1978) 
a 21K v2 


Bee 2 
+ doz * Ky Te Gee) see ee ial) 


YRR 
Cm saddeetaonal assitidption of neglecting the far field 
effects, 


< yea 


IYarl < o¥RI > 


2s the "parabolic" approximation for radial transmission 


[erocwl1I?e). The resuit is the Leontovich~Fock (1946) PE 





Ye = i1(AtB)y F 
where 
Ao == 
2K 327 
0 
and 











eech 2+ the form used by Brock (1978). 
be Implementation 

Brock (1978) used the SSFFI to solve the PE 
because it has several significant advantages which outweigh 
its disadvantages. Advantages include exponential accuracy 
=n depth, second-order accuracy in rang, energy conserva- 
tion, unconditional stability and computational efficiency. 
Disadvantages are a uniform mesh and periodic boundary con- 
Gditicns to satisfy the FFT, and filtering of the sound speed 
Meontle Dy Ssmocthing discontinuities <9 avoid spurisus high 
angular-frequency components (Brock, 1978). Additionally the 
algorithm assumes a flat pressure-rsiease surface, a vanish- 
ing field at the maximum depth and 32 pssuis-radiation condi- 
Som a= “he Peter-botto® interface by smoothly attenuating 
Eiemeecta (preck, 1973). 


A numerical algorithm cf 


vos i1CA+B) ; 
LARC Bt+E | 
WCRt+AR,Z) = BENE A) Ghee) 3 
leads =o 
7 AR = LAL - o 
Were = e882 erp yg at8** Pre (yCR,2)) } 
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OF 


LARK (n?-1)/2 ee IG 2 
W(RtAR,Z) =e FFT !{e OPPTCWCR,Z))} 


Moec@erlitate the computations since a spatial FFI and its 
snverse are required to transform between depth and wavenun- 
ber (Brock,1978). This procadure is implemented by two 
alternating steps, *+tke first of which considers propagation 
*m a WomoOgenectS medium +o account for iiffraction and the 
Beeena ce account for refraction (Brock,1978). In kis 
aGevelopment Brock gives an alternate calculation of the 
envelope function stemming from 


WCRAR.Z) = a tdRA/2 _tARB a tARA/2 NOR) 


Pals oyeeters slichtiy from th2 previous 2xpression contain- 


tJ 


Beg ¢he FFT and inverse Fr. hus i+ is ast currently known 
Namenmesetaticnship is used by Brock; this could be resolved 
Gagec by a detailed examination of the computer code. In 
any event as McDaniel (1975a) has shown, either approach 


would be of sufficient accuracy with appropriate selection 


Of AR and AZ. 
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ine Soup appears §*0 be quite satisfactory 
Pebetcically, but i= has to simplify the astual oceanic vari- 
ebility to allow realistic computation. It apparently will 
handle source and receiver depths and horizontal separation 
distances well as noted by its degrees of accuracy in depth 


and range. The smoothing of the sound speed gradients and 


Dek wom Srttreec=s will have to be considered. 


Gne LGeens Of chis descrip‘sioa has been to pro- 
vide a base from which comparisons can b2 made with other 
mime@isccal solutions and modifications +*5 this algoritha. 
The effect of the source model on the WI can be evaluated by 


+hese comparisons. The primary E€iements of the WT are the 


real and imaginary parts of th 


(p 


complex oressure amplitudes. 
Tiws the algerithmic differences with respect to computation 


@e che pressuee cield will be =he primary basis for the con- 


b ty 


Paltiscon cf effects on the WT. 
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McDaniel (1975a) compared th2 basic Tappert ana 
Meme Gin SSEPT,.:Which is used by Sroeck, *o9 normal mode theory. 
She found that errors arise in phase and group speeds, 
@itkough one mode can be propagated with the correct ampli- 
micencnad Group and phese spescis. McDaniel describes three 


sources of ezror: (1) approximzting the wave equation by th: 


SE, 





eet es Horec above, (2) limiting the range steps and per- 
mossible sound speed gradients, and (3) truncating the field 
Beue mate cepth. The result »9f the SSFFT will be a shift 
in the modal interference at long ranges if many modes are 
propagating (McDaniel,1975a). The error due to range step 
dimitations is third order in the range step increment 
(McDaniel, 1975a), thus a small range step is required to 
Mame Wace thes error. The third srror is related to inclusion 
of a highly absorbing botton which causes preferential 
eattenuaticn cr higher-order mnod2s; the dominant contribution 
Will thus result from the lower modes (McDaniel, 1975a). 
Brock @ al (1977) have described a procedure 
20 Zmprove this condition. It maps the index of refraction 
and depth such that «he phase Speeds and turning points of 
the nermail moc¢s are preserv2i. This involves constructicn 
of a “pseudoproblem" such that the parabolic phase speeds 
Well be eque. to the elliptic phase speeds of the corre- 
sponding modes in the original problem (Brock et al ,1977). 


This is given as 


, %, 


ec ee cu 


Coes 
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Tt is found to greatly improve the SSFFT agreement with the 
Sl iigtiic solution. This correction technique has one basic 
constraint, which is that isospeed regions are not allowed 
in the sound speed profile because of poor eigenfunction 
maoping (Brock et al, 1977). 


GQ. Subsequent Modifications 





ae Introduction 
There have been two significant modifications 

Since the originai algorithm was published by Brock (1978). 
These are source function and volume attenuation alteraticns 
(@@%20,1977) and the incorporation of a bottcem and variable 
Bamiee Seep (Stieglitz et al ,1979). The Tatro modification 
was developed in 1977 but was added to the Brock algorithno 
after 1978. Discussion of these two aspects will compiete 
mame @esersptacn of the PE algorithm installed at NPS and 
Meo tazed iw this study. 

Be.  Taero Modifications 


The, Gaussian Soutce originally utilized by Brock 


+e 


(1978) was consider)ed to be inefficient in terms of progran 
Sige, Ai@tfa-ion and an upper frequency limit which is zrela- 
tively Lom (Tatro, 1977). A new source function vas devel- 


oped. This was intended to solve the problem of large 


4 1 


saa ey 


tp 
yet 








vertical wavenumbers from a Gaussian sources that were previ- 
and isd to an 


range (Tatro, 1977) 
low pass filter in 
centered at the 


eliminated with 


ously 
emeG@rEpect depiction,of the field. 
is basically a 

is 


This source 
domain which 
initial vertical wavenumber 


*he verticai wavenumter 
Bewree depth (fatro,1977). Tha 
field will thus be a constant up to a prescribed value and 
zero ror higher wavenumbers. The \us2 "Of the filter iis 
yield improvement over the Saussian source by 
h2one POCccuS Waen teons orm ng 


claimed to 
Minimizing the aliasing which 
(Tatro,1977). The benefits derived from 
reduced computation time and 
frequencies 


<9 the depth domain 
the use of this new scurce ars 
Stor age and the availability fee gies 
(@atre,1977). 
He also “mOdified the Brock algorithm by intzro- 
ducing volume absorption as calculated from the equation of 
Gitore (1967). The decision can be made whether or not to 
inciude volume absorp<ion. 
@. Seicqlitz et almModi fications 
Broek (1978) had tonsidered the PE most useful 
>> lena targe prepagation at low frequencies ani smell 
grazing angles. The bottom was nodeled 2s fully absorbirg so 
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that the field would vanish at the maximum depth of the 


() 


transfcerm. This depth was obtained by extending the water 


column depth by one-fourth to obtain the transform depth. 


fhe inde) of refracticn in this region was 


where a and b are empirical constants. The energy in the 


bottom was attenuated such that any additional modes which 


j- 


resulted from transform truncation were minimized. This 
empirical expression was selected to coincide with 2 compa- 
rabie normal mode solution. This rather crude model of the 
bottcm appeared +o be acceptable for sitiaations wherein the 
bottcm hed a relatively small 2f fect. 

S2eegla- Zz etoel.(1979} considered it necessary 
Sov ette@pt +o develop a aore realistic ocean bottom since 
acoustic propagation was being utilized in more varied and 
extreme situations where botton interaction was important. 
They developed =wo additional options for handling <he bet- 
tOoa> (49 a pazrt#ally absorbing bottom with specification of 
reflective icss versus grazing angle of the equivalent ray 


and (2) the direct insertion of sediment scund speed and 
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attenuation profiles. The algorithm derives the sediment 
profiles from the grazing angle loss values in the firs¢+ 
opticn. The supvlied profiles are used directly in the sec- 
ema Option (Stieglitz et al ,1979). A new complex nodified 


index of refraction 


[Meee = {n-(R-Z) = 1} + 2p(Z)/K 
‘@) 


fs ceppaced if the bottom attenuation is to be calculated. 


ne ecadwtioral change, the calculation of the 


p 


Banmge st@p, has also been added by Stieglitz et al (1979). 
Tt consists of an initial search of the vertical wavenumber 


domain for the méximum pressur2 component. This is then fol- 


ts 
ct 


iswed by ase@rch for the firs componéen* which is 50 dB 


below the maximum component. The wavenumber of <his 


compcnent 


is svamed to require a range step of 


pls jelo=: ae, 
AR A/S 50 ' 


GU 





where 


oe 
0 fe) 


Tecwaange) Step was fixed in the current research for the 
Bemeenoe Of he spatial transform in the Wl, thus the above 
description was not utilized. It has been included to con- 
plete the description of the computer coisa which apparently 
comprises the NPS mod¢l, PEMODEL. 
5. SSEPFT Evaluations 

The SSFFT has been invastigated on several occasions 
tO @xamine specific items which might be modified to yield 
better agreement with actual propagation. Most of these 
Seeeeseaea) yeth the algorithmic description by Tappert and 
Ha@rdan (1973). Examination of these studies will note fur- 
*7WE>spossibie limitations on accuracy in the SSFFT which can 


xtended to tzhe aT. 


oO 
(D 
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Peemoereid (1975) examined the accuracy of the SSFFT 


as a function cf range and frequency. It was shown that the 


aj 


Safer G4 


(p 


Sonsc-ated greater accuracy at isnaer ranges, the 


foyer the rrequency. The SSFFT was considered a gcod 


approximation in =the examples given for 100 HZ and 10 HZ =o 





about 111 and 15000 km, respectively (Fitzgerald, 1975). 
The current study will examine transmission at 50 and 100 HZ 
at ranges to 100 nautical miles. These ranges will exceed 
the limits of Fitzgerald in some cases. 

McDaniel (1975b) investigated splitting the total 
field into transmitted and reflected fields in deriving the 
5 Sed oS It was necessary to neglect the reflected field for 
the matrix associated with the Tappert-Hardin SSFFT, since 
the two fields did not decouple wher the wavenumber was 
independent of range (McDaniel, Tab Thus an exact 
expression could not be determined for the transmitted field 
mecnowe whe reflected field. (The retlected field is neg- 
Meeeee =r tae Cirrent research since i+ uses the SSFFT.) 


This will iead to phase speed and group 


in 


peed errors as 
McDaniel (1975a) pointed out earlier, as discussed in S¢ec- 
Pivem LIt-p-c. These errors were Melamezea Dy Brock etal 
Se 7yeeey 2 Gdpoing technique for both tne index of tefrac- 
tom @@a ~iWe Gepth. ven so, the implication of this possi- 


© error cn the WF is that the interference patterns may 


p-* 


D 
pet "se Seee-, Which will lead to inexact or biased wavenua- 
oper nulls. McDaniel (1975b) orovosed another matrix which 


4 


achieved the Gesired decoupling and thus allowed expressions 


" 





or the transaitted and reflectad fields. This alternative 
w2S not investigated in this initial research. It may be a 
way to resolve non-exact results and could be studied later 
z0 isolate guantitatively the biases related to various 
input parameters for the WT. 

Palmer (1976) related possible inaccuracies in the 
SSFFT to appreximatiors involving relatively small terms in 
the exponential approximations for the opsrators A and B and 
also to the use of a single reference wavenumber. These 
appreximations will net be investigated specifically in this 
research but they may be related to potantial disagreement 
tecwmecn “ne SSPPT and tne theorstical Lloyd's Mirror effect. 
The associaticn of errors due to a singl2 wavenumber can be 
reccegnized =rom the wavenumber spectra, which will be pre- 
sented lecter in section IV-A. They will show not only an 
energy concentration at the ref2rence wavanumber but signif- 
*cant energy Gistributed, as well, over a range of values. 
Daimer (1976) emphasized that the physical properties of the 
ocean must be considered before terms can be neglectsd, such 
as the spreading of energy over several wavenumbers. 

DeSanto (1977) studied the relationship betwsen the 


pressure and the envelope function (called velocity 
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potential by DeSanto), which is directly relevant to the 
current research where pressure is calculated from the mnod- 


eled envelope function. DeSanto arrives it 


l-2a 


y(R,Z) = AR TZ.) Q(R,t,Z) exp (ik /Zt)(R2+t?) pee! oa 


where Q satisfies 


ae + 2{(i-a)/R *CIK RI/tIQ, + San + 2iK Q, + 2Q,3_{in(p(t,Z))} = 
9Z 


Ze 1 
K {K,(t,Z) ~ K,CR,Z) }Q 


sv 
{3 
Au 


Ko(R,Z) = KCR,Z) - Ky (Z)e_ 


ine 8h For =he*envelcpe function is obtained if this func- 
zion is ¢valuated by stationary phase; thus, it appears that 
a more accurate result could b2 obtained by better evalua- 
*ion of the integral. DeSanto 2t al (1978) developed a cor- 


ESGce tere ror =he SSFrtT, known as «he Corrected Parabolic 
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Approximation (CPA), that handled phase inaccuracies better 
than the SSFFT. Further study might address the response of 
the WE to the CPA with respect to improved phase accuracy as 
a function of range. . 
6. Environmental Sensitivity of the PE 
This final section in the analysis of the PE will 
consider the expected environmental sensitivity of the PE, 
and, by extension, the WT. The important relationship of the 
PE tc the WT is the input of the complex pressure amplitude 
and ultimately «he resultant wavenumber distribution. As 
discussed eariier in Section II-3, the wavenumber could be 
associated to the source angle and, by axtension, to the ray 
type or family (Lauer, 1979). Consequently, oceanic parame- 
Mees Whee woadity the transmission of c2rtain rays can be 
related to the resuiting wavenumber distribution. 
The sound speed profile will favor certain ray paths 
a= GetGgeoent depths ln the water column. There should be a 
corresponding shift in ray paths and hence wavenumbers as 
she source and receiver depth geometry is changed. Tae 


Geme@enana surface-refiected rays will travel at diferent 


speeds jiependin uponmche gmadiencts alomg the psmefile and 
hus lead to differen spatial interference patterns. This 


49 





can leai to wavenumber shifts and changes in the null spac- 


ties used in the current study will include real 


-~ 


ang. Prot 
examples and isospeed cases. 

Another important aspect is the transmission geon- 
etry. Generaily the model will be run at considerable range 
and in relatively deer water, which will be consistent with 
the applicability of the PE. A series of source depths with 
a single fixed teceiver depth will be run to generate ths 
wavenumber nuil spacing values that will formulate a deter- 


mination curve for each case of geometriztal variation. The 


j+- 


distripution of these resulting curves can be examined to 
evaluate the effects cf transmission geometry on the WT. 
Two frequencies will 02 examined in a sinilar man- 


ner. This would be one of the most important parameters to 


MWoeure sSstudts=s Since it @dulad have 32 direct 


tt) 


investigate in 
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(p 


aring on operational applications. It is expected that 


here will be differerces in the resulting wavenumber spec- 
Se2 “dependirg on the input frequency. mhe telatzonmski p 
between frequency and transmission geometry will aiso be 


studied since these are the easiest paraneters to vary dur- 


t’- 


1g a 


3 


2xperiment and will lead to the most complets anaiy- 


{nh 


2S CE @ocei Varaability. 
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The surface will be assumed to be planar and pres- 
Sure-release so that there wiil be complete specular reflec- 
tion. The minimum source depth used will be that at which a 
definable interference pattern san be observed for the spec- 
ified frequency. The slope approaches zero near the surface 


on the source depth determination curves. This result may 


}- 


fy 
i-! 


UGase wat the Wr in its current form is probably not 
applicable for sufficiently shallow sources. 

The final impcrtant area of variability will be the 
bottcm boundary conditions. Tfhere are several options now 
available in the PEMODEL. The fully absorbing and partially 
absorbing bottoms will be csonsidered. A <fcully absorbing 
bozttcm will pe used in most cases to simplify the interfer- 


ence pattern by elimirating any rays which interact with the 


bottom. Thus the basic interference pattern may be Evaiu- 
ated. Varicus bottoms can then be tested tc study the 
erregn’ ser che bottomon the dT. Differing bottom less 
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curves will be compared in ai <> demonstrate quali- 
Getewery thet «he bottom does have an 2ff2ct. The third bot- 
~om cpotiorn of @ realistic sedimant sound speed could also be 
investigated once the general effect of the bettom has been 


s 


@eteftmeneac. it woulé not offer any additional significant 
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results at this point beyond those realized from the fully- 
and partizaliy- absorbing bottoms. 

All input parameters ara related and thus the exami- 
mation of any one particular factor must be considered with 
respect to the full range of variation of the other factors. 
The number of combinations of parameters would be very large 
for a model such as the PEMODEL. Thus this research can 
only begin to explore general relationships. Additional 
research is necessary to delineate specific parameter 


oe 


SHpacts on the WT. 


C. FAST FIELD PROGRAM 





1. General Description 
The rast Field Program (FFP) was selected as the 
Meee Model <0 CoMpeare to the PEMODEL sincs it was the ori- 


G@2=n Gf the Wl and the initial testing of the NWT by Lauer 
(4779) Ses performed with the FFP. The description of the 
FFP which follows is summarized from DiNapoli (1971). 

Mae §“Dassac recult of the FFP is direct numerical 
eei@eieoesOn Dy the apflication 9£ the FFT *9 fieid theory to 
compute propagation credictions in a minimum computation 


time. The pressure cen be represented in the forn 
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where G, the Green's function, must Satisfy 


d2G + {K*(2Z) =k? 16 = Seep 
dz? 


and the associated boundary conditions. The field equation 


Can then be @rit+en as 


seme ee ae 
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Roeetess term, © Gan be obtained through use of the FFT at 
each renge incremert. Given p as a function of range, the WT 
eculd be applied as usual. DiNapoli compared an example 
Irom the FFP with the normal wmode progran of Bartberger and 
AcKkKler (1973) and found exceilenzt agreement between the 
locations of the peaks and the real parts of the eigenval- 
ues. Tim chiy significant difference is that the FFP 


includes more higher order noi2s but dots not include the 
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Tiest four sacdes. DiNapoli states that the FFP offers a 
Significant reduction in computation time yet still provides 
teasonable results as compared to normal node theory. This 
is accomplished by expressing the sound velocity profile as 
one or more exponential functions of depth. These exponen- 
tlai finctions allow the use of recurrence relations to 
quickly calculate the input to the FFT. 

The FFP offers two options for the modeling of the 
bottom; a two-layered fluid bettom with specified sound 
speeds and a semi-infinite bottom of one sound speed. The 
propagation effect Erom the first would be eventual return 


Geogr roam ths second bottom interface. The 


iD 


Of Lays emee refl 


semi-infinite bottom would not allow any return of sediment 
rays. 
¢- Comparison of FFP and 2E 


The FrP is a range indapendent addel and the PE is 
range dependent. The environnent was considered constant 
Geer Ghee renge in the current research, and thus the two 
models san be compared. The FFP uses the sound speed profils 
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erent approach 


tty 


Ss an Wepenential function which is a daéf 


the linear segmentation of the PE. The number of expo- 
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Sal ©unc=:ons could be sufficient t> model an actual 


54 








profile approximately as well as the PE. The PE bottom mod- 
eling does nct appear to offer any significant advantage in 
theory over the FFD, which allows for a semi-infinite or 
lavered botton. ieeis lNely £rom a thesretical standpoint 
ome tee PE and FFP inputs to the WE could be similar enough 
to allow comparison. The wavenuiber Spectrum shows a series 
of peaks which result from the nearby singularities of the 
Gaeen's function. These can be expressed as the normal 
Bodes. The utility of the WI stems from the spacing of the 


nulis between the envelopes which contain these peaks. 


De NORMAL MODE THEORY 
1. General Description 
The fundamental importance of ‘ths normal mode solu- 
taleme and the reason that is it used to check the validity 
Of oOtmer Models, is that it is an exact solution to the wave 
Sqma@tion (Rinsler et al ,1982,p.430). The following descrip- 


«ion continues from Kinsler e+ 21 (1982,pp.430-432). 


The solution for a point source is given by 


We 3+ / 6-92) p= 1 6(R-R.) Ss ) 
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Analysis yields 


- Jwt hae = | 
D e ( JT) RZ (20)2 (2)H (KR) 


and Zn must satisfy 
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dz? 


Maete Kh is constant and appropriate boundary conditions. 
PHos solutien is for trapped modes and does not include 
hose which are evanescent. Normal mode theory can be 
iinked to the WT by the spatial represantaticn of the wave- 


humber 
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The rays in the deep sound channel which have the same iocal 
G@-@6ction Or propagation correspond <9 the normal modes py 


“he reiatiorship 
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which extends to the normal mode phase speed 


where 1 designates the 1th mode. 

Normal mode theory lenis itself t9 many applications 
in underwater acoustics and it may be formulated for both 
range dependent and irdependent Situations. A dependent form 
which is appropriate for comparison with the SSFFT was 
developed by Kanabis (1975). He developed a normal-mode 
model which ailows large changes in the ssund speed profile, 


o. 
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Pthy 48a Dcttom composition with range. The range domain 
is segmented into regions within each of which sound speed 


amd bottom composition are functions or dtpth but not range. 


Q) 


The normel node solutions for #2ch region are matched at the 


(D 


2. Comparison of NM and P 


ju 


Te would be expected that the WE would respond simi- 
larly co these two modeis. This is based upon this example 
of @ normai mode orcgram and supported by the previously 


iscussed evaluation studies involving P& and NM. The NM 
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Bee on= ere generally appiicable to long range propagation 
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at low frequencies, which is quite similar +o the PE. The 
NM can be identified clesely with the WT ina physical 
sense, because of the interfarence effects in the NM that 
result from several multiple reflections (Officer, 1958,p. 
117). Officer further states that the rnuaber of reflections 


Moll Whereas and the time interveais between the incidence 


Or successive reflections will decrease as the range 
increases. This might be comparable to depth variations 
M2eh Will be shown for the WT. The NM ylelds exact solu- 


tions for modeling with planar boundaries. Thus the results 
of a WT based on NM would be limited only by approximations 
MeeGh Meee intrinsic to the WT. A comparison orf PE- and NM- 
based WP data should aid in evaluatin the oMpace, “Of tie 


Barcicubar PE algorit ha. 


2. PEre@ee DITPERENCS METHODS 
1. General Description 


One final group of methods which has evolved 


2, 


Fecencly ate the finite difference methods of solution for 
*he parabolic equations (Lee and Papad2kis, No 9). They 


EeeemBOre gemerai solution which will be appropriate in 


7) 
w 
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Shtiestlew water or wher¢ bottom interaction will be impozvan. 


Toey EWeener state that the FD and 9DE nethcecs We 124 © be 
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Cuemeretor <O SSFPT kecause an artificial bottom is not 
required for the benefit of tha FFT. Expl’ cat and Imp. text 


inite difference schemes are considered as well as an ordi- 


th 


nary differential equation. The explicit finite difference 


(EPD) scheme is 


where a” represents y(Rn,Z0), Ra is the nth range point, and 
7 2© the atk depth point. Dine’ anplteLt tliat e"ditierence 


(IFD) method from Lee and Botssas (1982) is 
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PWeEtED scheme “Gis consistent, stabie and convergent in the 
solution of the PE (Lee and Papadakis, 1979). 


The ODE method is 


where a and b are related by 


/ 
fond 


wee = a(K ,R,2)u + DCK RR, Zu, : 
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Mies ctecpnigue is stable, consistent ard convergent also 


Which Suppor=s its u 


F)) 


e. On2 additional aspect of the ODE 
approach is the calculation of 2 variable range step. 
¢- Comparison of PD and ODE with PE 

Omen fieouit of the FD and ODE “asSthods is greater 
fiexibility in arbitrary bottom and surf2zce boundary condi- 
tions than that available with the SSFFT. Lee and Papadakis 
(W979) State that while the FD and ODE methods are both 
superior to the SSFFT, the bes= appears t> ba the ODE, fol- 
towed by the EFD, which is limited by ier Thus the FD 
and ODE methcads appear to be significant improvements over 
the SSFFT. Tkis result should extend to 2 WT based upon FD 
or ODE input. 

More recentiy Lee et al (1981) have re-examined the 
IFD and rourd it to be as equaily efficient as the ODE. The 
generai advanzage cf the IFD over the SSFFT is that the 
problem is solved within the water columa since the botion 
eztects can be modeled by straight line segments. This is 
SeQ@eter =o the ODE, but the IFD doés not require =he amount 
Of storage necessary for «he current ODE (Lee st al, 1981). 
McDaniel and Lee (1982) have extended the IFD to the treat- 


Water ot Werz~.:cai density discontinuities. They found zhat 
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while without this more realistic treatment of the water- 


ox 


Sitom interrace the IFD approximates the SSFFT and differs 


n 


Mag@2crcantly from normal mode, the IFD corresponds 
extremely well with the normal mode calculations with the 
interface. In conclusion, th2 numerical tachniques related 
to FD and ODE appear to be viable alternatives to «he SSFFT. 
They will permit a more accurate representation of underwa- 
ter acoustic transmission, which will be directly relatable 


to the WT. 
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A. GERERAL DESCRIPTION 

Twenty scenarios were selected to examine <=he general 
aspects of variations in several oceanic and geometric 
paramet2rs. A set of input parameters anda sound speed 
profile were constant in each scenario while the source 
depth was varied. Five to ten source deoths were used for 
each situaticn. The vertical and horizontal axes of the 
wavenumber spectra correspond *> the spectral intensity and 
“ke scaled wavenumber beta, r2spectively. The intensi-y 
axis Wes normelized to unity for ease in plotting. The null 
spacing was determined for é€ach spectrum. These null spac- 
“ng values were plotted versus source depths =o yield the 
determination curve icr ¢ach scenario. 

Table 2 iists the parametric specifications for each 
SCE Reario. Bogece 1G «he can be divided into isospeed or 
sound speed profile cases. An isospeed profile of 4890 ft/s 


he 
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ScpOe@ PRs Swsfece to the botton Was ceiec'ad *o evaluate 


ometry and the general consistency of experi- 
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Bonmeae cesui-s, assuming a fully-absorbing bottom. The aT 





= based On the @erivals of the direct and surface-reflected 
rays; ‘thus, the isospeed cases will be the easiest to conm- 
pare with theory. There will be ro sound speed differences 
as a function of depth which could modify the arrival times. 
Additionally, several realistic ocean profiles were studied 
+o allow a qualitative comparison with theory and explore 
spacial and *emporal variability. It is recognized that the 
theoretical Lioyd's Mirror 2ffect applies only «+o isospeed 
profiles. Results from real profiles were compared with this 
theory to see if they could approximate it. The spatial and 
<empceral scales of variability used in this research were 
approxinately 200 miles and 20 days. 

Mors ané@éiyses were performed wit ‘he isospeed cases 
Bane@e thos was an initial invastigation. It still remains 
to consider other simple cas?s, such as linear negative and 
positive gradients, surface ducts and de¢p sound channels. 
Ris@ realistic ocean studies that would involve greater *en- 
poral and spatial variations would be beneficial. Time was 
Migited diirirg =his research, thus efforts were divided 
BetWeen BoGei validatior ani 2ceanic sensitivity. i ees 


béli¢ved, however, that ths ccmparisons which will be 


(p 
t-- 


pue@ieeeea are sufficient for a preliminary investiga*tion 


SiO themcenS tivity ci the WT. 
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TABLE I 


Scenario Speci fications 


ESCSTEED PROFILE 


Run Rov Depth Range Frag Water Depth Number of 


Seok BDgett R onm Ponz Gees SD Runs 
HO 300 50 109 1000) 7 
r1 300 50 50 10000 6 
T2 300 100 59 10000 6 
pe 300 50 100 5009 6 
T4 300 100 100 5000 5 
I5 300 50 50 5009 6 
T6 300 100 50 10009 6 
E7 300 100 100 10000 5 
r8 800 50 50 10009 5 
jigs) 10000 50 50 10009 5 
R1 300 25 109 10009 5 
R2 300 50 100 10009 

R3 300 75 199 10009) 5 


REALISTIC PROFILE 


Run Rov Depth Range Ea a a ti SEG PSern 


Set Rn Dae T KR 0m us lle SD(number) 
AQ 300 100 1090 9000 8 
Al 1000 100 100 9009 sy 
AZ 9000 100 190 9009 5 
A3 300 100 100 9009 5 
BO 300 100 10 12009 8 
CO 300 100 100 9009 8 
DO 300 100 109 12009 8 
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The compiex pressure input data were not modified so 
+hat the entire wavenumber spectrum would be observed at the 
F2ceiver. These erfects would need to be evaluated if the 
WI was to be compared to received signals. The surface and 
bottom were assumed flat. Boundary variations would have to 
be considered for WT use under realistic ocean conditions. 
A nominal beam widthof thirty degrees was used and in 
almost all cases the range step was 0.1 mils. 

Frequencies of 50 and100 Hz, water column depths of 
5000 and 10,000 feet and ranges of 50 and 100 nautical miles 
were the basic comparisons for most isosoeed scenarios. A 
receiver depth of 300 feet and a fully absorbing bottom were 


M@eidized, wniess otherwise noted. Additional comparisons 


= 


ere performed, such as range variations to 25 and 75 miles 


fi) 
J 


G@ ueecsimer depth *ariations of 800 and 10,000 feet. A 
memage Of 100 mautical emiles and a Erequency of 100 Hz were 
need for =he actual cases unless othérwis= noted. Th bottom 


d 


wv 


Dube were 9000 feet for profiles A and C and 12,000 feet 


EOr @erot—2f#es B and D. 


Gige peeiivstic profiles wers obtained @rom the Acoustic 
Storm Transfer and Response Experiment (ASTREX) which was 


Gonew@esea by SPS Qn the northeast Pacific Ocean in November 





and December, 1980 (Dunlap,1982) . The ASTREX study was con- 
ducted to evaiuate the effect of oceanic response to atmos- 
pRette StTOLTSS Or acoUstic propagation. The profiles were 
selected for this research <>) represent only spatial and 


temporal variability without respect to any experiment oper- 


sv 


e2.0nS’. The profiles were obtained from expendable bathy- 
thermographs (AXBT) which were dropped by U. S. Navy P=-3 


e@trecrett. Fig. & Shows the relative locations and times for 


ct 


these profiles. The digitized profiles were analyzed by the 
University cf Hawaii to remcev2 spikes and false values. 


TWeSe profiles were combined with climatological salinit 
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ry 
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Pectiwe= at the Fleet Numerical Oceanography Center in Mon- 


Tabie II gives the d¢scriptive data for éach 
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prozile. The bottom loss curv>s which were used are similer 


j 


To) tome £OUnd in Urick (1979) and were arbitrarily devel- 


Opedeeo am@zcate higher loss for profiles 8 and D. This is 
Befteree To HbCttoOm cofposition; bottom typ2s for profiles A/c 
and B/D aze calcareous sand and Clay; respectively 
(NA VOCEANO,1978) . Thus nigher loss is expected for B/D. 


Meral feature of the spectra for almost all test 


ty 


Aq 
Gtoee ts a ad¥stinct U-shape wir peak values at the minimun 


and maximum wavenumbers. Figs. 5 through 10 show this phe- 





Pigucte 4. 


Space and 
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relat 
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TABLE II 


Lower en Prtosmatten for Realistic Profiles 


Peo file Tua t Long Date Time (Z) MLD(m) 

A G1-24N 127-23HW 16NOV 1928:04 SZ 

B YS-47N 133-05W 16NOV 2121:25 61 

G 41-23N 128-22W obEC 1ehSs.37 58 

D US-46N 133-098 oseeEC 2083345 79 
romenon for isospeed cases and Figs. 11 through 16 fora 
pediieqe.¢ seurd speed profile. The abbreviations on these 
figures are source depth (SD), receiver depth (RD), and 


wa zer column dep*h (Z) in feet; range (R) in nautical miles; 
and frequency (F) in Hz. This key will also apply to the 
source depth determination curves which will be presented 
Waser. Tke ~igh- and left intansity maxina appear to corre- 
Spond td beem eleveation anglzs cf 0 and 30 deg, respec- 
sel y, @Whech could be considered to approximate the direct 
and surrace-~rerflected waves. A 30 deg beam angle was used; 
Rous, tnere Way be a relationship with the left maximun. 
This beam angie is also near the maximum angle which Tappert 
(W677) cohstGereq to he appropriate. Ther the left maximun 
Bight be a result of algorithmic inaccuracy. Different bean 


eng@im@s ccoul@ be tested £0 deztarmine if this would explain 


os 





the unexpected U-shape curve. This apparently anomalous 
spectral patzern might be explained by a possible transforn 
assumption of symmetry about zaro and th2 inclusicn of neg- 
ative transtorm values. itwes also possible that the left 
wavenumber peek may be the result cf aliasing of the high 
wavenumber values. DiNapoli (1971) suggested that aliasing 
Goma GcCcur 2nd thus he utilized only one-half of the ranges 


Meaela. nas has not been tested in the cu 


4 
4 


ent research; 


however, examination of different range séeaqments might 


resolve this apparent problem. A remot2 odsssibilizy is that 
the direct ard surface-reflectsd rays each exhibit 2 Bessel 


functies d2splays This might be depictsi by a peak at th 


iD 


rh 
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reference wavenumber for each ray with a set of interferin 
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envelopes becween the peaks. A mathematical explanation for 


(D 
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Phis @Osscibility has not yet been pLOr 
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ay No spikes are 


ce 


(D 


S2en 42h Wery sheliow Water as depicted in Fig. 17, but also 
mo Gigemact “mul pattern is pbserved, thus near surface 
source depths have been neglected. 

These figures were also included t2 show a series of 
source depth runs and tne resulting decr2ise in beta spacing 


as a fumetien of seurce depth. The null spacing was manu- 


ally measured for each run as the distanc= between two adja- 
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RUN 143 
9D:200FT RD:3O00FT F:100HZ Z:SOO0FT R: 
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| NORMALIZED SPECTRAL INTENSITY 
0.000 0.025 0.050 0.075 O.100 QO.325 O.150 0.1758 0.200 
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SCALED WAVENUMBER BETA (1/FT) 
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RUN 144 


S0:300FT RO: 300FT F:100HZ Z:SO00FT R: 100NM 


NORMALIZED SPECTRAL INTENSITY 
0.000.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 
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RUN I45 
SD:SOQOFT RD:3O00FT F:100HZ Z:SO00FT R:100NM 
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RUN 146 
SD:800FT RO:SOOFT F:100HZ 2:SQ00FT R21 





NORMALIZED SPECTRAL INTENSITY 
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SCALED WAVENUMBER BETA (1/FT) 


Figure 8. : Speed @zorile Ravenumber Spectrum for SD of 
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RUN 147 
SD:i1S500FT RO:300FT F:100HZ Z:SOOOFT R:100NM 


NORMALIZED SPECTRAL INTENSITY 
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RUN 148 
OFT RD:3O0OFT F:100HZ Z:SO00FT ayia 
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RUN AQS 
SO0:200FT RD: SOOFT F:100HZ Z:SQ000FT R:100NM 
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RUN AQ4 
SD: 300FT RD:3OQOFT F:100HZ Z:90Q00FT R:100NM 
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RUN AQS 
SD:SOOFT RD:300FT F:100HZ Z:SQ00FT R:100NM 


0.4 0.5 0.6 0.7 0.6 


0.3 


, if 


{dil ll Wy add | \ 


0.050 0.048 0.040 0.035 0.030 0.0200 0.035 0.010 0.005 6.000 
” SCALED WAVENUMBER BETA (1/FT) 


: 


NORMALIZED SPECTRAL INTENSITY 


0.2 


0.0 0Q.!1 





eee Sees tic Profile Wavenunber Spectrum for SD of 
| Le 


718 





ral 
} 6. 


Q 


M 


RUN AQ6 
SD:800FT RD:3O00FT F:100HZ Z:SOQ00FT R:100NM 
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RUN AQ7 
SO: 1SQOFT RD:300FT F:100HZ Z:9000FT R:100NM 
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RUN 14 I 


SoU: SOFT RO: SO0FT Fs: 1Q0HZ Z:SOQGOFT R:1CQNM 
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SCALED WAVENUMBER BETA (1/FT) 
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Gomeeieeeno1 -y iminima, such as shown in Figs. 5 through 16. 
The nulls selected for measurement wert those in «*h= middle 
Se Ghee Spectraey ewhich led t5 a consistent technique for 
comparison and avoided the rapidly decreasing null spacing 
near the left intensity maxima. The distance was measured 
=> te neazes= 1764 of an inch and equatei to the respective 
mull spacing. Thus, there is some inheren*+ measuremer= 
SanOmem ete wexec- null becations were difficult tc datermine 
in the shallew isospeed cases because of gent+le spectrun 
CagWetmee ana upcertainty in auli location. The oe 
prefiles were very ccmplex and again the nulls were diffi- 
cuit to leeate. One value for ¢ach specz=rum plot was usad 
=> geneeate a point on the determination curve, such as Fid. 


13) Ecce a2 See eos of cms. Th] variation »f the data points 


j— 


ino mw 2=ee cheorsctica 


a 


Mowe SweaarOSeecusve Willa pe <=he normal 
distance between zhe foint and the curve. Measurement inac- 
curecy will decrease as sourc2 depth increases. Th ss wea 


be siognificar= pecause at shallow depth there is some varia- 


tion noced from sound speed differences. The mui iS?) fats 
Wel db-@ere neo 22 areet depths, but there is little variation 


already, thus the measurement errors were net considered 


Significant. It is possibie that a numerical algorithm which 
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Gols “Guenti~atavely identify the null locations could maini- 
mize this errcr. 

This expectation is inferred from the wavenumber plots 
by Lauer (1979) and his examination of range determination. 


He rtound that range determination resolution increased as 


the distance between source and receiver decreased. hens 
pOssibde thet a2 SHRiiar relationship might exist for depth 


déterminatior. Thus greater resolution nay be realized at 
shallow source depth. Even at jreater source depths the res~ 
clution may be acceptabie for certain applications where an 
approximate depth is acceptabie. 

Teeseeed'c-ion curves which result from «he individual 
Mr Specera wiil be discussed in foilowing sections. Table 
Sil =s% iicting of the prediction curve comparisons which 
Mere conducted. In the figuras only on? symbol may appear 
2&2 sO@eee depth :f beth curves had the same null spacing 


value. 
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— SOURCE DEPTH DETERMINATION CURVE 


SOURCE DEPTH (FT) 
500.0 1000.0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 


0.0 





0.0000 0.0025 0.0050 


0.0075 0.0100 0.0125 0.0150 90.0175 0.0200 
BETA NULL SPACING (1/FT) 


Pogue elon) sxa@ple of Source Depth Detsrnination Curve. 
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TABLE Lf IT 


Comparisen Specifications 


mS OSEESO PROFILE 


Figure Run Sets Constant Parameters Varying Parameter 


=e Sa ee ee ee ee ee ee ee a ee ee ee ee a Ss Dw ee 


13 Pit > Weare... F 250,RD: 300 Z:10000,5000 ft 
20 Dt 8 goa, ©: 100,RD 2300 Z:10000,5000 ft 
21 i2,7uo KE 100,F:50,R8D 2300 Z:10000,5000 ft 
ge Tine Mm: 100, F>100, RB: 300 2:10000,5000 ft 
23 >, mer F > 5092 :5000, RD: 300 R3:50,100 nro 
24 ys PF so, a 10000, RDS300 Reo OperOO nt 


Peewee eee) 6 100, 27170000,R2D; 300 R:25,50,75 am 


aif I2,I7  R:100,Z:10000,RD:300 F:100,50 hz 
27 ii 1 22°50 sae YOO00, RD:300 Eero =O iz 
28 Tdeat 8,19 R50) £250.25 1.0000 RD: 300,800, 10000 ft 


ReALTSTEC PROFILE 


Figure Run Sets Constant Parameters Varying Parameter 


Ze nO WBO) R: 100, F:100, RD:300 Space:200 nao 

30 Soyo Rs TOO Pere PRD: 300 Space:200 nan 

31 eorCcO R: 100, F2 100, RD: 300 T=nres20°> days 

ee £0,D0 3: 100,F:100,RD: 300 Time:20 days 

33 80,43 Fs 100, 16NO¥, RD: 300 R:100,50 am 

36 <AOgA1,A2 7F: 100, 16NOV, R:100 2D:300,1000,9000 f= 
COMPOSITE 
Biguewe® pun Sets Constant Parameters Varying Parameter 


35 (Ali Tsospeed and Actual Profiie Run Sets) 
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Be is0Ss PEED CASES 
1. Watez Column Lept! 
Mog. Voets tet CO@perisen for a frequency of 50 Hz 
and a range of 50 n@ at water column depths of 5000 and 
10,000 feet. It was not expected that tnere would be any 
VaTiation if the water cclumn depth was modified, since all 
bottcm interecting rays should be fully absorbed. The 


curves showed excellent agreement at source depths of 800 


feet and greater, but there was some variation noted for 500 
f2et and shaliower. Some variation may be related to meas- 
urement error, as noted abov2. Maximum variation from 


“Heomr “es observed for both depths in the middie depth 


Tang€, approxtmeteiy 500 feet, Which is O= 


ey) 
}-s 
W” 
O 
ct 
2 
m 
4 
(D 
Q 
1 
O 
an | 


greatest curvature in the sourc® depth determination curve. 
Fig. 20 is the same 2xcenpt the frequency has been 
inere@sed to 100 Hz. There is better agreement from 200 to 


3000 feet; nowever, the significance of she frequency varia- 


raj 


(4) 
rh 


tion has not been evaluated. ge 2 1 or a frequency of 
50 Ge apaea serge of 100 crap, 2nd Fig. 22 zs for 100 Hz and 
100 nm. Some variaticns are noted; however, they may not be 


Sige Wetoeme, Wai-ch would agrees with the Supposition that 
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SOURCE DEPTH DETERMINATION CURVE 
TSOSPEED WATER DEPTH VARIATION 
| RsSONM FsSOHZ RD: 300FT 


2900.0 


LEGEND 
© = LLOTOS THIRRCR 
@- Ji: Z=-10000FT 
O- 15: Z-S000FT 


SOURCE DEPTH (FT) 
1000.0 1500.0 2000.0 

e 

* 


$00.0 





0.0 


0.0000 0.0025 0.0050 0.0075 0.0100 0.0125 9.0150 9.0175 4.0200 
BETA NULL SPACING (1/FT) 


Sow re 19. re ees Profiie SD Determination Curve: Z 
ees et ton — 1. 


— al 
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differences will not be significant for 


water column depth 


.sospeed profile. 
deerecuency sof 50 HZ 


100 no. 
The 


fe 
IS 


an 


2. Remge Vesdati c 
23 1s the comparison for 
ROG 50 and 


of 5000 feet 
between the curves. 
3000 feet but 


water column depth 
=o 
from inaccura- 


and a 
There is essentially no difference 
experimental curves agree well from 800 
once again tne shallow variations may result 
om the spectrun plots. Pig... cous 
deoth has been 


water column 
increases the frequency 


cies in evaluating beta 
comparison, but the 
feer. Hig. Vio 

50 and 75> files. Again all 


which indicates 


a Similar 

increased to 10,000 
tom Cam@es or 25 
With theory, 

SOUreCS wall prob 


£0 100 Hz 

curves show very good agreemen* 
eyond 100 nm were not tested 

SSFFT was not a 


Ranges b 
the 


Tange variations within 10) nm of «the 
ed thax 
This could be 


thet 
abiy be insignificant. 
Since Fitzgerald (1975) suggest 
Coeamcpomnotcmme. On £Or 100 HZ be yond “11 lekao. 

of 50 nm and 


apa 
a range 


Llawer. 
a 


tested 
HZ. 


3 Hi 
® > -_ 
——— 


equen 
-“s the colfiparison for 
or 50 and 100 


feet, 


deptn of 5000 feet 
well with theory below 800 
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SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED WATER DEPTH VARIATION 
R:SONM F:100HZ RD: 300FT 


7300.0 3000.0 
: 


LEGEND 
@ = LLOYDS MIRROR 
w- 10: 2-10000FT 
aie. a= J3: Z-SOOOFT 


| SOURCE DEPTH (FT) 
1000.0 1500.0 2000.0 


300.0 





0.0 


0.0000 0.00205 0.0050 0.00775 0.0100 0.0125 9.0150 0.0175 0.0200 
BETA NULL SPACING (1/FT) 


m= gure 20. ESosPe rd Proz=te SD Detarmination Curve: 2 


Sriaw on - A 
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SOURCE DEPTH (FT) 
2000.0 2500.0 3000.0 3500.0 4000.0 


1300.0 


Figure 21. T 
3 V 


SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED WATER DEPTH VARIATION 
R:1Q00NM FsSOHZ RD: SOOFT 


LEGEND 
@ = LLOYDS MIRROR 
s- 12: Z-10000FT 
a o- 16: Z-SOOOFT 


0.0000 0.002 0.0050 0.0075 0.0100 0.0125 0.0150 0.0175 
BETA NULL SPACING (1/FT) 


= 
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bation 


G1 





0.0200 


cspeed BOSSE SD Determination Curve: Z 
b Gm = 3s 





| SOURCE DEPTH (FT) 
500.0 1000.0 1500.0 2000.0 2500.0 3000.0 


SOURCE DEPTH DETERMINATION CURVE 
_TISOSPEED WATER DEPTH VARIATION 
R:1Q0NM F:10Q0HZ RO: 300FT 


LEGEND 
® = LLOYDS MIRROR 
e- 17: Z=10000FT 
be o = 14: Z=SOOOFT 





9.0000 0.002 0.0050 0.0075 0.0300 0.0125 0.9150 0.0375 0.0200 


BETA NULL SPACING (1/FT) 


Prof-le SD Seas trms2racion Curve: Z 
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= SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED RANGE VARIATION 
FsSOHZ Z:SO000FT RD:3S00FT 


2300.0 3000.0 3500.0 4000.0 
e 


LEGEND 
e- LLOYDS MIRROR 
wp- 15: R=-SONM 
Q- 16: R-100NM 


SOURCE DEPTH (FT) 
1500.0 2000.0 

e 
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1000.0 


900.0 





0.0 


0.0000 0.0025 0.0050 0.0100 0.0125 0.0150 0.0175 0.0200 
BETA "NULL SPACING (1/FT) 


Figure 23. isospeed Profile SD Determination Curve: & 
Vamiatszon - 1. 
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SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED RANGE VARIATION 
FsSOHZ Z:10000FT RD:300FT 


2300.0 3000.0 ¥500.0 4000.0 
e 


LEGEND 
e - LLOYOS MIRROR 
we ji: R=SONM 
Ge ]2: R=}QONM 


SOURCE DEPTH (FT) 
1000.0 1500.0 2000.0 
= 


500.0 





0.0 


0.0000 0.002 0. 


0.00775 0.0100 0.0125 0.0150 9.0175 90.0200 
BETA NULL SPACING (1/FT) 


rigure 24. Isospeed Frofile SD Determination Curve: R 
ha Sect2On w+ 2. 
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SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED RANGE VARIRTION 
F:1Q00HZ Z:1Q0Q000FT RD:300FT 


3000.0 
© 


2900.0 


LEGEND 
e- LLOYDS MIRROR 
@- Ri: R=-2SNM 
O = R2: R=SONM 
o = R3: R=-7SNM 


SOURCE DEPTH (FT) 
2000.0 


1500.0 
a 
« 





0.0000 0.0025 0.0050 0.0075 0.0100 0.0125 0.0150 0.0175 0.0200 
BETA NULL SPACING (1/FT) 


Peom@e@e 2s. Esospeed Frofile SD Determination Curve: R 
Peeler yon - 3. 
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feet the higher frequency appears to agres better. Fig. 27 
extends the water column depth to 10,000 feet and there is 
epparently better agreement with the theoretical curve. 
Overall there appears to be better agreenent at 100 Hz than 
ex 5® Hz, but this result ls not certain with the limited 
nuaber Of cases run. 


4. Receiver Depth Variat: 


ti: 
IO 
i 


Fig. 28 compares receiver depths of 300, 800 and 
V0, ©0Gmeecet for a range of 50 nm, water column depth of 
10,000 feet, and frequency of 50 Hz. Again all curves agré¢e 
reasonably well. There is perhaps slightly better agreement 
for the shaliow receiver, but t+his suppositicn cannot be 


Vet ifie@ Erom these Limited dati. 


5. Analysis of Isospeed Ca 


tw) 


es 

Thepe do not appear to be any significant differ~ 
ences for water column depth, range, frequency or receiver 
Geren @Fer <<he ligzts of variations considered. This was 


expected since the theory is supposed 9 be only a function 


~s 


: thus che we would 


(D 


O 
thy 


Se@mece depth for ar isospeed profil 


expected to be equally effactive for varicus transmis- 


LO 
(v 


(iD 


sion geometry patterns and frequencies, iz the sound spee 


oa 


ig2 is constant and the bottom fully absorptive. There 


rth 


oO 


22 
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es SOURCE DEPTH DETERMINATION CURVE 
TSOSPEED FREQUENCY VARIATION 
R:100NM Z:10000FT RO:300FT 


3000.0 3500.0 4000.0 
° 


2500.0 


G- 17: FeiQ0HZ 


SOURCE DEPTH (FT) 
2000.0 


LEGEND 
e@- Loe MIRROR 
2: F=SOHZ 


1000.0 1$00.0 
© 
& 


$00.0 


0.0 





0.0000 0.0025 0.0050 0.0075 0.0100 0.0125 0.0150 0.0175 9.0200 
BETA NULL SPACING (1/FT) 


Eemeate 26-5 beosoneed Profile SD Determination Curve: F 
Variation -1. 


ed! 





4000.0 


SOURCE DEPTH (FT) 


SOURCE DEPTH DETERMINATION CURVE 
TSOSPEED FREQUENCY VARIATION 
~ ReSONM = Z:10000FT RD: 300FT 


2500.0 3000.0 3500.0 
6 


= 
: LEGEND 
@- LLOYOS MIRROR 
r=) e- ]}: F-SOHZ 
: a a= 10: F-100HZ 


1000 .G 


$00.6 





0.0 


0.0000 0.0025 0.0050 0.0075 0.01/00 0.0125 G.0!18SO 0.0175 0.0200 
BETA NULL SPACING (1/FT) 


Zi. “Seeepeea Freofiile SD Dezermination Curve: F 
Variation -2. 
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SOURCE DEPTH DETERMINATION CURVE 
ISOSPEED RCV DEPTH VARIATION 
R:SONM F:SOHZ 2Z:10000FT; 


2900.0 3000.0 300.0 4000.0 
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LEGENO 
e - LLOYDS MIRROR 
@- 1): RO-300FT 
Q~=- 18: RO=-800FT 
° = 19: RO=-1000FT 


SOURCE DEPTH (FT) 
1000.0 1s00.0 2000.0 
: 


| 
| 


300.0 





0.0000 0.0025 0.0050 0.0075 0100 «= 0.0125 (0.0150 «0.0175 = 0.0200 
BETA NULL SPACING C1/FT) 


Po gure 28. 1 eee Peeei te Sp) PEczermM:naeticen Curve: RD 
ari aye 


me ate ee le 


a9 





may be slight effects for receiver depth and frequency, but 


these nave not been substantiated. 


S. RPoALISTIC PROFILE CASES 
Page c? ise ths COMparison Of Spatial variability for 
meme CURVe= which resulted from the spectra of profiles A and 


B. (The sound speed profiles corresponding to the pre- 


bee 


diction curves will be included for each comparison.) These 
were both sampled on November 16, 1980, and represent a spa- 
tial separation of abcut 200 nm. Both profiles result in a 
Ste2ght bas below tne theoretical Lioyd's Mirror curve. 


Ttece se More Warjetion eiong each curve, which 1s to be 


expected because the profiles are not isospeéed. Also, par- 


{~4 


Saaeee Y ea SOLT D2 h bettoms are now used 2s compared to the 


uily abserbing bottoms used with the isospeed casés. The 


rh 


tf) 


rey) 


1£¢ 


co 


'ty 


nd bottom effects hav2 not yet been separated to 


fu 


etermine which contributes more to the departure from 


theory. This would require an 2xtensive s2t of intercompari- 


(n 
O 
oe) 
(Nn 
On 
@ 
ct 
2: 
10) 
@m 
;3 
2 
3 
D 
Nn 
D 


shese two parameters. (Wes "D5ESeENt Lntene 2s 


j-4- 


ko ShoweeWat cc@an variability is importan:.) Pieees, aes 
cult to difierentiate between the two curves, but it appears 


Beomewewe 25 a Sltqhe trend in 38 ‘tsward smaller null 


100 





Spacings. Sound speed profile B is slower than A over the 
entire depth. A slower sound speed would appear to indicate 
that a ray of lesser angle would be required to provide the 
second arrival for the interference pattern. That would cor- 
respend to @ greater vertical wavenumber and decreased beta 
null spacing. 

Fig. 30 compares profiles C and 0D, which were san- 
pled at the same locations as profiles A and B, respec- 
tively, but on December 6, 1980. A similar pattern is noted 


With D spacing less t~ 


id 


lan C, which may again correspond to 
Slower sound speeds fcr profils D. 
2. Temporai Variation 
Fig. 31 compares oprofiles A and C which represent 
“n2 Same lecation, tut aseparation of 20 days. There 
Somearewes be 4 Curve crosSing at about 8)90 fee~, such that 
Speceamg Ss Wess at A for shallow depths and greater for A 


pelcw 800 feet. The A sound spp2ed profile is slower than the 


C profile down to 800 feet and the pattern is reversed frcm 
ou ZO B2000 zee. Slowex speeds appear +o be related to 


* 


snotter muli svacings, whic 


a 
jp 
Ae) 


Mea=sS, cDaSiStent With “his 


example. The variaticns between the curves, however, are 
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SOURCE DEPTH DETERMINATION CURVE 
REALISTIC SPATIAL VARIATION 
Fs1QQOHZ R:100NM RO: 300FT 
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= QO- PROFILE 8 
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Vania icom @8,2): SD Determination and 
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SOURCE DEPTH DETERMINATION CURVE 
REALISTIC SPATIAL VARIATION 
F:iQ0HZ R:100NM RO:300FT 


4000.0 


LEGEND 
@- LLOYDS MIRROR 
@- PROFILE C 
o- PROFILE D 
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2000.0 3000.0 
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1000.0 
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BETR NULL SPACING (1/FT) 
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Figure 305 Spatial Variation (C,D): SD Determination and 
AXBT Curves. 
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SOURCE DEPTH DETERMINATION CURVE 
REALISTIC TEMPORAL VARIATION 
Fs1QQHZ R:100NM ROD: 300FT 
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1 Variation (&,C): SD Determination and 
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Fig. 32 presents a similar comparison for profiles B 
and D. The speed profiles show almost no differences, except 
eRa:e D is Siightiy slower. This result appears as a 
Slightly sharper deep sound channel which will mean slightly 
ello ver welocity f05 “the direct path. Thaz effect may be 
reflected as a slightly greater null sgpacing as shown. 

Results for profile A were compared for rangés of 50 
and 100 Sh@® in Fig. 33. NO significant variations were 


noted, which is Similar to the results for *he isospeed 


cases. 
4. Rece-ver Depth Variation 
Restiits for profile aA were also compared for 


Eecemye: dep-hs of 300, WOOO and, 2000 “fest an Fig. 34. 


There may be slightly better agreement for the shallow 


(p 


Cem@tee, bus “this is not cerzain. 


D. GORPOSITE COMPARISON 

Fig. 35 is a composite of all isospeed and actual 
profiie data sets. Tsospec¢i data show slightly better 
agreement wich *«heery; however reaiistic profiles ars a rea- 
sonakl¢ approximation to the theoretical curve. Realistic 


profile variations from the isospeed curve appear to be 





Sadecea tc the scund speed profile and different bottom 
boundary conditions. The general relationships of the date 


indicat2 thac below about 800 f2e* ro significant variations 


are noted. Specific curves will apparently provide better 
predictions at shallower depths. I+ may be sufficient fer 
Geet caublscat10ns <O provid>= climatological or theorsti- 


ed 


c21 predicticn curves. These curves could be enhanced of 


modified by updated curves that would reflect the current 


sound speed profiles. 
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SOURCE DEPTH DETERMINATION CURVE 
REALISTIC TEMPORAL VARIATION 
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Pigute 32. Zemporal Variation (8,D): SD Determination and 
AABT Curves. 
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SOURCE DEPTH DETERMINATION CURVE 
REALISTIC RANGE VARIATION 
F:1Q00HZ RD: 300FT 
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LEGEND 
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Piguasess3. Range (Rk) Variation for Profile A. 





SOURCE DEPTH DETERMINATION CURVE 
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Figure 34. Receiver Depth (RD) Variation for Profile A. 
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SOURCE DEPTH DETERMINATION CURVE 
overeeo ANO REALISTIC CASES 
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Fogure 35. Ccmposite Comparison for Isospzed and Realistic 
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The Wavenumber Technique (WT), as proposed by Lauer 
(1979) has been investigated for source depth determination. 
Several idealized and actual scenarios of oceanic and acous- 
+ic parametric variations were tested. There is excellent 
agreement between the ideal test cases and the «theoretical 
Lloyd’s Mirror effect. The raalistic sases showed signifi- 
cant similarity as well, but with some bias toward smaller 
wavenumber spacing. This is qualitatively attributed to 
sound speed profile variations in th2 water and sediment. 

There appears to be a response of the WT to variations 
in 2mpetjy wnich is beneficial in allowinzy 2valuation of the 
relative daportance of individual parameters. Changes in 
range and water column depth do not indicate any Significant 
Yariation. Some differences are possibly suggested for 
receiver depth and frequency, but the results have not been 

° 
fully analyzed. The WT appears *o9 respond well +o changes 
“nN the EeRvironmental parame=3rs, suco as sound speed 


pao £ ¥ic. This is seen from th? average bias cf the actual 


(D 


ar +o be 


cases from «he isospeed cases. The bias does 


po 


PP 





relatively consistent and shows the importance of oceanic 
variations in both space and in time. Thus, it appears that, 
ur inteor@Mation is available for the bottom structure, and if 
@arealistic sound speed profil? can be generated, then a 
very specific family cf curves could be generated for vari- 
ous frequency and receiver depth combinations. This simplic- 
ity indicates possible GeeLaetenal applzcabizlity and 
flezrwpai it 7. 

These are the first known source depth determination 
curves tor the WT and thus ther2 is no conparison available. 
It may be that there are some biases in these curves which 
would b2 discovered when comparative runs are2 made using 


other inputs ocr voropagation nod2ls. Comparison with observed 


ct 


Signals might also reveal biases tha ace related to the WT 
or to ‘he prepagation nodel. An explanation has not been 


determined for «he U-shape nature of the wavenumber spectrum 


plots and thus slightly different curves might result if 


SE 


<haG is found to bias the nuli sp2cing. Also there is a 


[ese Weerod GF Marualiy estimating the 


ct 


bias associated wi 
pull spaceng. ThiS @ignt be r2solved with 2 computer proce- 
dumee” =O Locate the au ls. The experimental curves actually 


ae eweWaen checry under the asSuhoption of a flat, 
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sotally-reflecting surface anda fully~-absorbing bottom. 
The premise that the theory might be approximately accurate 
in realistic cases appears valid. This overall generally 


Emoeth £it suipperts the validity of the WI. 


i-- 


This has been an initial investigation of the WI, and 
thus cnly a broadbrush study was possibls. These prelimi- 
nary results indicate that additional analysis might be ben- 
eficial. The response of the WT to input from different 
propagation models would deternine if results from different 
models could be compared. Detailed investigavtion of the spe- 
CaGic etfects of varicus varameters would lead to a quanti- 
SearcarWe atalysis cf the sensitivity of the WT. An explanation 
tor the amomalous U-shaped svectrum and the development of 
an algorithm which can accurately determine null locations 
would be necessary before the WT could be used practically. 
An extremely important investigation would be the cemparison 
of analyzed actual signals +t9 predictions from the WT =o 
GEtermine applicability, Secrecy, consistency and any 
apprcpriate biases. 

The intert of this research was to datermine =he valid- 
-zy of the WT by development of source iepth determination 


curves and to evaluate the oc2anic and acoustic 





sensitivities of «hese curves. The results obtained indicate 
that in general the WE appears to be valid wi*h respect <«o 
=né bhboyd's Mirror effect. The parametric variations which 
were investigated show that oceanic variability is extremely 
important, but that the WT is flexible for other parameters, 
Such as transmission geometry. The WT aopears to be viable 
for eperetionel applications. The additional investigations 
suggested above would determine the valiiity of this Suppo- 
Sition. It is conceivable that the WI could be developed 


nee @n cpeneational prediction product. 
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